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Abstract 

We present a calculation of the magnetic moments for the baryon octet and decuplet 
using the background-field method and standard Wilson gauge and fermion actions 
in the quenched approximation of lattice QCD. Progressively smaller static magnetic 
fields are introduced on a 24 lattice at beta=6.0 and the pion mass is probed down 
to about 500 MeV. Magnetic moments are extracted from the linear response of the 
masses to the background field. 
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1 Introduction 



Magnetic moments are an important fundamental property of particles. They 
determine dynamical response of a bound system to a soft external stimulus, 
and provide valuable insight into internal strong interaction structure. Ef- 
forts to compute the magnetic moment on the lattice can be divided into 
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two categories. One is from form factors which involves three-point func- 
tions [1,2,3,4,5,6,7]. The other is the background field method using only two- 
point functions [8,9,10,11]. The form factor method requires an extrapolation 
to Gm(Q 2 = 0) that is potentially fallible due to the discrete momentum on 
the lattice [12]. The background field method, on the other hand, is free from 
this problem since a static magnetic field is applied to the gauge background. 
Here we report a calculation in this method. It is an extension of our work on 
electric and magnetic polarizabilities [13,14,15]. 

2 Method 

In order to place a magnetic field on the lattice, we construct an analogy to 
the continuum case. The fermion action is modified by the minimal coupling 
prescription 



where q is the charge of the fermion field and is the vector potential de- 
scribing the background field. On the lattice, the prescription amounts to a 
modified link variable 



Choosing A y = Bx, a constant magnetic field B can be introduced in the 
z-direction. Then the phase factor is in the y-links 



Dp = + gGf, + qA i 



(1) 



(2) 



Uy B ^ = exp (iqa 2 Bx). 



(3) 
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On a finite lattice with periodic boundary conditions, to get a constant mag- 
netic field, B has to be quantized by the condition 

?5a 2 = ^, n= 1,2,3, ••• (4) 



to ensure that the magnetic phase factor is periodic in the x-direction. How- 
ever, for N x — N y — 16 and 1/a — 2GeV, the lowest field would give the 
proton a mass shift of about 390 MeV, which is too large (the proton is 
severely distorted). In this work, we abandon the quantization condition and 
choose to work with smaller fields. To minimize the boundary effects, we work 
with fixed (or Dirichlet) b.c. in the x-direction and large N x , so that quarks 
originating in the middle of the lattice have little chance of propagating to the 
edge. 

We use 24 4 lattice at f3 — 6.0, and six kappa values k=0.1515, 0.1525, 0.1535, 
0.1540, 0.1545, 0.1555, corresponding to pion mass of 1015, 908, 794, 732, 
667, 522 MeV. The critical value of kappa is k c =0.1571. The strange quark 
mass is set at k=0.1535. The source location for the quark propagators is 
(x,y,z,t) = (12,l,l,2). We analyzed 87 configurations. The following five dimen- 
sionless numbers r] = qB a 2 =+0.00036, -0.00072, +0.00144, -0.00288, +0.00576 
give four small B fields (two positive, two negative) at eBa 2 =-0. 00108, +0.00216, 
-0.00432, +0.00864 for both u and d (or s) quarks. Note that "small" here is 
in the sense that the mass shift is only a fraction of the proton mass: \iBjm ~ 
1 to 5% at the smallest pion mass. In absolute terms, the field is enormous: 
B ~ 10 13 tesla. 
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The mass shift in the presence of small fields can be expanded as a polynomial 
in B, 

5m(B) = m(B) - m(0) = c ± B + c 2 B 2 + c 3 B 3 + c 4 B 4 + ■■■ (5) 



To eliminate the contamination from the even-power terms, we calculate mass 
shifts both in the field B and its reverse —B for each value of B, then take the 
difference and divide by 2. So in terms of cost, our calculation is equivalent 
to 11 separate spectrum calculations: 5 original r\ values, 5 reversed, plus the 
one at zero field to set the baseline. In light of future dynamical simulations, 
the factor can be reduced to 3 if only one nonzero field is desired. Another 
benefit of repeating the calculation with the field reversed is that by taking 
the average of Sm(B) and 5m(—B) in the same dataset, one can eliminate the 
odd-powered terms in the mass shift. The coefficient of the leading quadratic 
term is directly related to the so-called magnetic dipole polarizability (f3). 

For a Dirac particle of spin s in uniform fields, 



where the upper sign means spin up and the lower sign means spin-down, and 
\i = g^s. In our data, m + and m_ correspond to the '11' and '22' diagonal 
components of the baryon correlation function in the absence of the back- 
ground field. They are the same within statistical fluctuations. The average of 
the two components is usually used to extract the proton mass. We tried the 
following three methods to extract the g factors, and found they are equivalent 
within statistical errors. 



E± = m± ± \iB 



(6) 



9 = 



( 



m + + m 



2m + m 



) 



(E + — m + ) — (E_ — m_) 



eBs 



(7) 
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m + (E + — m + ) — m_{E_ — m_) 
9 ~ eB~s 

(Ej - m%) - (Ej - mj) 
9 2eBs 



The results quoted are from Eq. (7). 



(9) 



3 Results and discussion 

Fig. 1 displays a typical effective mass plot at the six kappa values for the 
strongest magnetic field for the proton. Good plateaus exist for all six quark 
masses. Our results are extracted from the time window 12 to 15. Fig. 2 
displays a typical mass shift, defined as 5 = g(eBs) from Eq. (7), as a function 
of the field for the proton. There is good linear behavior going through the 
origin for all the fields when the quark mass is heavy, an indication that 
contamination from the higher terms has been effectively eliminated. This is 
also confirmed numerically by the smallness of the coefficient in the B 3 term 
as shown in the same figure. At the lightest quark mass (lower right corner), 
there is a slight deviation from linear behavior. For this reason, we only use 
the two smallest field values to do the linear fit at all the quark masses. 

Fig. 3 shows the results for the proton and the neutron as a function of pion 
mass squared 1 . Note that the g factors directly extracted from the data are 
in the particles natural magnetons. To convert them to the commonly-used 
nuclear magnetons (hn), we have scaled the results by the factor 938/M where 
M is the mass of the particle measured in the same calculation at each quark 
1 Since pion mass squared is proportional to quark mass in QCD (m^ oc m q ), it is 
equivalent to plotting as a function of the quark mass. 
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mass. The line is a simple chiral fit using the ansatz 
jj, = do + aim w . 



(10) 



There is a lot of effort in the literature on magnetic moments and their chi- 
ral extrapolations based on effective field theories [17,18,19,20,21,22,23,24]. 
Although different approaches produce the same chiral behavior near the 
chiral limit, there are issues surrounding the range of validity and model- 
independence away from the chiral limit. The lowest pion mass in our results 
(around 500 MeV) is probably too large for a meaningful application of the 
results from these studies. Our goal here is to present the lattice data. The 
simple ansatz in Eq. (10) serves only to show that there is onset of non-analytic 
behavior as pion mass is lowered, so a linear extrapolation is probably not a 
good idea. To get some idea about the systematic uncertainty in the chiral fit, 
we tried two other different forms, one is 

\x = a + aim^ + a 2 m^ (11) 

and the other the Pade form [24] 

/i = a /(l + aim n + a 2 ml). (12) 

The extrapolated value to the physical point is 3.04(6) for the proton and 
-1.84(3) for the neutron from Eq. (10), 3.78(5) and -2.03(2) from Eq. (11), 
and 4.21(5) and -1.84(2) from Eq. (12). 

Fig. 4 shows the results in the octet sigma channel. The perfect agreement 
for the £~ should be taken as a coincidence. Fig. 5 shows the results for the 
four delta states. The agreement with the well-known Q~ moment in Fig. 5 
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is particularly encouraging since no adjustment is needed for this point (it is 
a direct comparison with the experimental value). We take it as a sign of the 
correctness of our calculation. The experimental numbers in the delta channel 
are not well-established. The experimental value for A + (slightly shifted for 
better view) is taken from Ref. [25]. For reader's convenience, all of our results 
are summarized in Table 1. Finally, Fig. 6 shows the proton and A + together. 
The opposite curvatures are a signature of quenched chiral physics. A similar 
behavior has been observed by [26,27] using the form factor method. 

About the issue of finite-volume effects as related to the Dirichlet boundary 
conditions, we did some further study. We repeated the entire calculation on a 
smaller lattice 16 3 x 24, while keeping all the paraemeters the same as on the 
larger lattice 24 3 x 24. So the only difference is the spatial box size (16 3 vs. 
24 3 ). Assuming the lattice spacing is the same (a=0.1 fm) on the two lattices, 
the dimentionless combination m n L x for the lightest pion is about 6.2 on the 
24 3 x 24, and 4.2 on the 16 3 x 24. The effects on the magnetic moments 
of the proton and neutron are shown in Fig. 7. Relatively significant finite- 
volume effects are revealed, even at the heavier pion masses, 'spoiling' the 
trend towards the physical point. We take this as indication that the 16 3 x 24 
lattice is too small. To really pin down the finite-volume effects, we should 
repeat the calculation on a larger lattice, say 32 3 x 32. But it is currently 
beyond our computing resources. 

4 Conclusion 

In conclusion, we have computed the magnetic moments of the baryon octet 
and decuplet on the lattice using the background field method and standard 
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lattice technology. Overall, our results are consistent with experimental ob- 
servations. Detailed comparison with experiment must await a full account 
of systematic errors present in the results, such as finite-volume effects. In 
addition, there is a need to push the calculations to smaller pion masses so 
that reliable chiral extrapolations can be applied. Nonetheless, our results 
demonstrate that the method is robust and relatively cheap, as compared to 
the form factor method. Only mass shifts are required. This may facilitate the 
push to smaller pion masses, perhaps with the help of chiral fermions (overlap, 
domain- wall, twisted mass, ...). Finally, we foresee no technical problems in 
doing dynamical (full-QCD) background-field calculations in order to remove 
the effects of the quenched approximation. 

This work is supported in part by U.S. Department of Energy under grant 
DE-FG02-95ER40907, and by NSF grant 0070836. The computing resources 
at NERSC and JLab have been used. 
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Fig. 1. Effective mass plot for proton mass shifts as the strongest magnetic field in 

lattice units. Plotted is the combination {E + — m + ) — (E- — m_) in Eq. (7) at all 

six quark masses in the order of decreasing quark mass from top to bottom. 
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Table 1 



The computed magnetic moments for the baryon octet and decuplet in nuclear 
magnetons (/tjv) as a function of the pion mass. The extrapolated values are based 
on Eq (10). The experimental values are taken from the PDG [16]. 



n 


0.1515 


0.1525 


0.1535 


0.1540 


0.1545 


0.1555 


Extrap. 


Expt. 


m w (GeV) 


1.015 


0.908 


0.794 


0.732 


0.667 


0.522 






p 


1.48(3) 


1.63(4) 


1.80(6) 


1.90(8) 


2.02(10) 


2.34(17) 


3.04(6) 


2.79 


n 


-0.94(2) 


-1.02(2) 


-1.13(3) 


-1.18(4) 


-1.25(5) 


-1.39(8) 


-1.84(3) 


-1.91 


£+ 


1.57(4) 


1.70(5) 


1.85(6) 


1.93(7) 


2.02(8) 


2.23(11) 


2.87(3) 


2.45 




U. OZ(l ) 


U.OO^l J 


U.Ot)(Z ) 


U.Oi-(Z ) 


U.Do^Z J 


U.D 1 \o) 


U. IO{L ) 


U.DD 


£- 


-0.54(4) 


-0.60(4) 


-0.68(6) 


-0.73(6) 


-0.78(7) 


-0.92(9) 


-1.48(5) 


-1.16 


E° 


-1.07(4) 


-1.10(4) 


-1.13(4) 


-1.15(4) 


-1.17(5) 


-1.21(7) 


-1.37(1) 


-1.25 




-0.76(6) 


-0.77(7) 


-0.77(8) 


-0.77(8) 


-0.77(9) 


-0.78(11) 


-0.82(1) 


-0.65 


A 


-0.59(2) 


-0.60(2) 


-0.60(2) 


-0.61(2) 


-0.61(2) 


-0.62(3) 


-0.70(1) 


-0.61 


A++ 


3.35(7) 


3.68(9) 


4.06(12) 


4.28(14) 


4.50(16) 


4.92(28) 


5.24(18) 


4.52(1.00) 


A+ 


1.55(4) 


1.67(6) 


1.77(8) 


1.80(10) 


1.80(11) 


1.64(24) 


0.97(8) 




A 


-0.002(0) 


-0.003(0) 


-0.004(1) 


-0.005(1) 


-0.007(1) 


-0.011(6) 


-0.035(2) 




A~ 


-1.58(4) 


-1.73(5) 


-1.89(7) 


-1.98(8) 


-2.07(10) 


-2.34(17) 


-2.98(19) 




£*+ 


1.40(5) 


1.56(6) 


1.72(8) 


1.80(10) 


1.86(11) 


1.88(17) 


1.27(6) 




s*° 


-0.13(1) 


-0.09(1) 


-0.03(0) 


-0.01(0) 


0.03(0) 


0.08(2) 


0.33(5) 






-1.68(5) 


-1.76(7) 


-1.83(9) 


-1.87(10) 


-1.89(11) 


-1.93(16) 


-1.88(4) 




E*° 


-0.09(2) 


-0.06(1) 


-0.(^2(5) 


-0.002(1) 


0.018(6) 


0.05(3) 


0.16(4) 
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Fig. 2. Proton mass shifts as a function of magnetic field B in lattice units at the 
six quark masses (lightest in the lower right corner). The slope of the mass shift at 
each quark mass gives the g factor corresponding to that quark mass. The line is a 
fit using only the two smallest B values. 
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Fig. 3. Magnetic moments (in nuclear magnetons) for the proton and neutron as a 
function of m^. The 3 lines are chiral fits according to Eq. (10) (dashed), Eq. (11) 
(dotted), and Eq. (12) (solid). The experimental values, taken from the PDG [16], 
are indicated by the empty symbols. 
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Fig. 4. Magnetic moments for the octet E. 
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Fig. 5. Magnetic moments for the A states. 
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Fig. 6. Comparison of magnetic moments of the proton and A + . 
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Fig. 7. Finite-volume effects on the magnetic moments of the proton and neutron. 
The solid (red) symbols are from the 24 3 x 24 lattice, and the empty symbols are 
from the 16 3 x 24 lattice. The experimental values are denoted by the stars. 
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